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Abstract 



The collective motion of a finite nuclear system is investigated by numerical simulation and by linear response theory. 
Using a pseudo-particle simulation technique we analyze the giant resonances with a multipole decomposition scheme. 
We examine the energy and the damping of different giant collective modes and obtain the dependence of these 
quantities on the proton-neutron ratio. In simulations of finite nuclei including only mean fields the centroid energy 
of the resonance decreases with higher asymmetry due to the change of the compressibility and the damping increases. 
The collisional correlation in turn leads to a decreasing damping with increasing asymmetry. 

Alternatively the giant collective modes in asymmetric nuclear matter are investigated within linear response theory 
including the collisional correlations via a dynamic relaxation time approximation. For a multicomponent system we 
derive a coupled dispersion relation and show that two sources of coupling appear: (i) a coupling of isoscalar and 
isovector modes due to the action of different mean-fiolds and (ii) an explicit new coupling in asymmetric matter due 
to collisional interaction. We show that the latter is responsible for a new mode arising besides isovector and isoscalar 
modes. The comparison with simulation results as well as experiments is performed. It is stressed that the surface 
effects and the collisional correlations are both essential in order to describe the correct damping behavior. A model 
is presented which allows to combine surface and collisional effects. Higher-order modes like the isoscalar dipole mode 
which where recently measured are discussed within this frame. 

Collective; motion beyond the linear regime is demonstrated for the example of large-amplitude isoscalar giant 
octupole excitations in finite nuclear systems. Depending on the initial conditions we observe either clear octupole 
modes or over-damped octupole modes which decay immediately into quadrupole ones. This dependence on initial 
correlations represents a behaviour beyond linear response. 
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Chapter 1 

Introduction 



Giant resonances are currently enjoying a resurgence of attention as a tool for the investigation of many particle effects 
in finite quantum systems. The study of giant resonances offers the possibility to learn more about nuclear forces. Of 
special interest are modern experiments with exotic nuclei which broaden our knowledge in a new degree of freedom: 
the variability in the proton to neutron ratio. The treatment of collective modes in nuclear matter is documented in 
an enormous literature starting from their discovery up to recent discussions |^ -[^5 and citations therein. 

We want to investigate the collective excitations in asymmetric nuclear matter |^ by numerical simulation of 
the Vlasov equation and by comparison with the linear response theory. While the former leads to an insight into 
finite-size effects like surface, the latter allows to consider collisional correlations in a control-able way. Both methods 
are contrasted and we demonstrate that the collisional correlations as well as the surface effects are important to 
describe the experimental damping of giant resonances. 

The damping mechanisms of collective motions in excited nuclei are a topic of continuing debate. Mainly two lines 
of thought are pursued. In one line of thought it is assumed that collisions are the only physical reason for damping 
which is described via a Fermi liquid approach (2^, The other line of thought considers new features of 

the finite nucleus, such as surface oscillations and a level density with finite spacing. The investigations are performed 
without inertia ||2l[- ||2^, ^ or by including inertia [^,^^- js^; note that inertia is absent in infinite matter. 

Both classes of models predict a comparable degree of damping necessary to reproduce the experimental data. 
Consequently, it is an open question which is the correct physical reason for damping. Of course, the correct description 
has to assume a finite nucleus consisting of nuclcons which are bound via the mean field, through which the nuclcons 
undergo mutual collisions and where the surface is formed by the particles themselves. These features are principally 



included in Boltzmann-Uehling-Uhlenbeck-(BUU) simulations ^ |4C[| or in its nonlocal extensions |41, |4^. In 



full simulations, however, we will not gain a simple insight into the physical origins of the damping mechanism, in 
particular, how much is contributed by the surface and how much by collisions. 

One aim of this article is therefore to compare both pictures in the frame of linear response theory. Within the 
collision-free Vlasov equation the linear response of finite systems is well known and allows one to calculate the 
strength function of finite nuclei. The resulting damping, however, does not reproduce the experimental damping of 
giant resonances since collisions are absent. This motivates us to develop a linear response theory including collisions. 

While most of the theoretical treatments of oscillations rely on the linear response method or RPA methods, large 
amplitude oscillations require methods beyond this level. In particular the question of the appearance of chaos has 
recently been investigated Q, ^ The hypothesis was established that the octupole mode is over-damped due 



to negative curved surface and consequent additional chaotic damping |4^, 48 , ^ . Here we want to discuss in which 
conditions one might observe octupole modes at least in Vlasov - simulations of giant resonances which will turn out 
to be dependent on initial conditions and are consequently an effect beyond linear response. 

The outline is as follows: In the second chapter we will give the numerical results of pseudo-particle simulation of 
the Vlasov equation. This will provide us with some insight into the magnitude of finite size effects on the damping. 
In the third chapter we present the linear response result of the Vlasov equation including collisional correlations. 
We will consider a simplified picture of infinite matter response and will consider the Steinwedel Jensen picture. We 
discuss two possibilities to include surface contributions in the linear response formalism. This allows us to describe 
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the experimental temperature dependence of the damping of giant dipole resonances as well as the structure function 

of the isovcctor dipolo resonances. The surface consideration allows to describe the isoscalar dipole resonance as a 
higher order mode. The comparison between linear response results and simulations is performed and finally we discuss 
nonlinear effects beyond linear response for the example of giant octupole modes. 



3 



Chapter 2 

Pseudoparticle simulation 



We will describe the giant resonance first by a kinetic equation using a pseudo-particle simulation ||30| . The kinetic 
equation for the quasiclassical distribution function reads for neutrons (for protons analogously) 



/„(p, r, t) + — 9r/n(P, r, t) - dr{Un + C/oxt)9p/„(p, r, t) = Icon- 

m 



(2.1) 



with the coUisional integral /corr discussed later and the self-consistent mean-field potential U given by a schematic 
Skyrme type |^ 

t/„/,(p,/) = a{f^)+b{fj±cl{f^) (2.2) 

with a — — 356MeV, h — 303MeV , c = 54MeV and s — 7/6. Here the neutron excess is / = e,^-Qp ^j^^ neutrons feel 
the H — sign potential, the protons the opposite one. 

By multiplying the kinetic equation by l,p or i? = ^ -f U respectively one obtains the balance for particle density 
p, momentum density u and energy density £ . Since the collision integrals vanish for density and momentum balance 
we get the usual balance equations 



dt dvj 

where the mean field energy of the system varies as 



dp{r,t) , d f d^dE 

i (2^ aF^^P' = ° 



dt drj {2Tihf dp 
dui{r,t) d f dp dE 



(27r;i)3 ,5/(p,r,t) 



^(:^: + W(P,r,t) (2.4) 



such that from (2^) follows the total energy density as 



{2TTn) 

With the help of this quantity the balance of energy density reads from (|2.1| 



d£(rA) d f dp ^dE , 9 ^ , , ,„ 

^ ' :E—f{p,r,t)^--E,orr{r,t) (2.6) 



at drJ {2TThf dp ' ' dt 
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with the correlation energy |£^, |5^ arising from the colhsional side. The coUisional side we wiU consider later in 
linear response. Here for the simulation we will neglect the collisions and will restrict ourselves to the Vlasov kinetic 
equation. In this way we will learn what are the effects of finite size and what are the effects of collisions. 



Let us now first look at the input for the simulation. The mass number dependence of binding energy £ /A from (2.5) 
is shown in figure 2.1. One sees that with increasing asymmetry the binding energy becomes weaker and consequently 



the compressibility decreases. The compression modulus is defined as the derivative of the Energy per particle 
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Figure 2.1: Left: Mass num- 
ber dependence of the energy 
for different neutron excess / = 
eji-eiL_ The dots denote the 

Qn+Qp 

saturation densities. With in- 
creasing asymmetry the satura- 
tion density and the binding en- 
ergy decrease. Right: Asym- 
metry dependence of compress- 
ibility modulus. With increas- 
ing asymmetry this modulus de- 
creases. 



,d{£IA) 



K^^n^^ (2.7) 



and took for J = the value K = 200 MeV according to (2.5). Figure |2.l| shows the asymmetry dependence of K. 
One sees that for neutron matter no binding energy would occur and a very low compressibility. 

We solve now the kinetic equation by representing the distribution function /(p, r, t) by a sum of N pseudo-particle 
distributions 

AN 

/(p,r,t)«/o(p,r,i) = ^_/s(p-pi(t),r-ri(t)) (2.8) 

i=l 

and use Gaussian pseudo-particles 

/s(P - Pi, I- - i-i) = c e-(P"Pi)'/2^? ^-(r~r^? /■2<yl (2.9) 
at ri with momentum pi ||5^. These pseudo-particles follow classical Hamilton equations 

p, = -VC/, r, = ^ . (2.10) 
m 

The only fluctuations introduced are due to some small unavoidable numerical noise |5^. We are using 300 pseudo- 
particles per nucleon and a pseudo-particle width of Ur = 0.466 fm. They are adjusted both in such a way that 
the experimental energy of the ^"^Ca giant monopole mode is reproduced. With these fixed two parameters the 
experimental behavior of centroid energy with mass number is than reproduced over the full range for giant monopole 
and giant dipole resonances. We have checked different numbers of test particles. The dependence of observables on 



the width is discussed in 1 56 . Numerically the ground states of nuclei are realized by Wood- Saxon shapes of density 



and Fermi spheres in momentum. 



2.1 Multipole analysis 

Provided we have now solved the kinetic equation we will have the distribution function represented by N pseudo- 
particles. Since we are interested in moments of the distribution function, the density, current and energy, we can 
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expand these moments in the test particle representation as well. The momentum distribution of these moments can 
be obtained by spatial integration: Fa{p) — /dra/(p,r) with, for the mass distribution, a = 1, for isospin, a — r, for 



kinetic energy, a 
they read 



2m 



, and for kinetic isospin energy, a 



2m 



respectively. Decomposed into spherical coordinates 



AN 

E 

1=1 



Radial integration determines now a spherical distribution 

AN 



1 Sy't 



(2.11) 



(2.12) 



which can be decomposed mto spherical harmonics 



oo I 



1=0 m=-l 
AN 

air. 



AN ^ 



(2.13) 
(2.14) 



i=l 



The observable distributions (p) are normalized to \/ An aoo, i-C. to mass number A, total isospin T, kinetic 
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Figure 2.2: The picture shows for isoscalar 
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monopole oscillations of IpZr the time evo- 



lution of multipole moments ^J^^^aiQ{t),l = 
0,1,2,3 corresponding to monopole, dipole, 
quadrupole and octupole oscillation for the 
distribution of isospin a — r, kinetic isospin 

2 

energy a = mass a — 1 and kinetic en- 



ergy a 



p 

2m ■ 



energy Ekm and kinetic isospin energy EkinT , respectively. The polar angular distribution of moments reads now from 
(|t|) and ( pT3| ) 



AN 



1 Si^~A) 



-N sin(i?,) Jo 



2ir 



dip Paid, ip) 



(2.15) 



;=o 
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with Legendre polynomials Pi. As a measure for the strength of the resonances we obtain now the coefficients aio of 
the corresponding moment from (2.15) as 



21 + 1 21 



^ AN 



An 



]_ Pi {cos d,) 
^ N°" sin(?9,) ■ 



(2.16) 



are displayed as a function of time in figures 2.2.2.5,2.8 and 2.9 



These amplitudes of multipole moments, a/o w 4^ 
The value oio means the dipole moment, vanishing for isoscalar resonances, 020 is characterizing the quadrupole 
oscillations and 030 the octupole ones, etc. 



2.2 Isoscalar monopoles 



The first analyzed mode is the isoscalar giant monopole mode which plays an important role in the determination 
of nuclear compressibility. The connection of the compressibility and the energy of giant monopole resonances is 
discussed e.g. in [57|. 

The figure 2.2 shows simulation results for a monopole oscillation. Here the excitation has been performed by 
adding an extra momentum to the test-particles in the direction of the center of mass. We have excited in this way 
a clear monopole breathing mode which can be seen in the oscillation of the kinetic energy. The corresponding mean 
field energy performs the opposite oscillations that the total energy is constant. The fact that energy is oscillating 
between kinetic and correlational ones describes why we have here a compressional or breathing mode. All other 
modes remain unexcited. The finite value of the isospin and mass for I = reflects the conservation of isospin and 
particle number. 



With the previously chosen fixed width of test-particles we see from figure 2.3 that the experimental mass number 
dependence of monopole oscillation is well reproduced over the whole range of mass numbers. However, the experi- 
mental damping can be seen to be largely underestimated by the Vlasov simulation which yields about 2 MeV. This is 
a first indication that the mean field, even for a finite system, cannot account for the whole damping and dissipation. 
Instead we have to take into account coUisional correlations which will be performed later. 
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Figure 2.3: The mass number 
dependence of giant monopole 
oscillation energy and damping 
width. The solid line denotes the 
experimental values while the 
dashed line describes the Vlasov 
simulation results. 



It is now instructive to examine the isospin dependence of the isoscalar monopole. Since this is the exper imen tal 
value which determines the nuclear compressibility, the isospin dependence is of direct importance. In figure 2A the 
dependence on asymmetry of the isoscalar monopole energy is plotted for a hard as well as a stiff equation of state. 
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Of course, the hard equation of state leads to a higher monopole energy corresponding to a higher compressibihty. 
Analogously to figure the monopole energy decreases with the asymmetry. The hard equation of state leads to a 
more linear decrease while the soft equation of state remains almost unchanged up to a certain asymmetry and then 
decreases faster. 
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Figure 2.4: The giant monopole energy (a) 
and normalized to the symmetric case (b) ver- 
sus asymmetry for two different param eteri - 
zations of the mean field according to ( ^.6| ). 
The stiff equation of state (black line) leads 
to higher monopole energy and shows a more 
linear decrease while the soft equation of state 
(grey line) shows a weak dependence for small 
asymmetry and decreases faster for higher 
asymmetries. 
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2.3 Isovector dipole resonances 



On the next figure |2.5| the multipole analysis for ^qZt oscillations is performed. The excitation is now chosen as an 
isovector dipole one created by a shift of proton against neutron spheres. One recognizes that a clear dipole aio mode 
is excited which is seen in the oscillation of the isospin and isospin energy. The kinetic energy now remains constant in 
contrast to the isoscalar mode because we have no compression mode. Like the monopole case we now investigate the 
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Figure 2.5: The time evolu- 
tion of multipole moments for 
the isovector dipole excitation of 



40 analogously to figure 2.2. 
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mass number dependence. With the same parameterization as for the monopoles, the experimental mass dependence 
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of the energy is well reproduced, figure 2.6. The damping is again under-predicted. The Vlasov simulation leads to 
a nearly constant width of 2 MeV. This is again a hint that collisions cannot be neglected. The isospin dependency 
of the dipole mode is plotted in figure 2.7 and is characterized by an increasing width and decreasing centroid energy 
with increasing asymmetry. 
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Figure 2.6: The experimental mass 
number dependence of isovector giant 
dipole oscillations (upper curve) is com- 
pared with the Vlasov simulation re- 
sults. The oscillations in the experi- 
mental damping width are due to shell 
closures. 
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Figure 2.7: The isospin dependency of the centroid energy and the damping of the isovector dipole (left), isoscalar 
quadrupole (middle) and isovector quadrupole mode (right). 



2.4 Quadrupole resonances 



Let us now investigate the quadrupole resonances. For these we distinguish isoscalar and isovector modes. The 
quadrupole resonances are excited by dividing the spatial distribution into two equal pieces which are accelerated in 
opposite directions. In the first example the neutrons and protons are in phase giving rise to an isoscalar mode. 

We see from the kinetic energy in figure |2.S| that besides the clear isoscalar quadrupole mode there is also a weak 
isoscalar monopole mode excited. The isospin dependency of the isoscalar quadrupole mode is plotted in figure 2.7. 
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As in the case of the isovector dipole mode one sees a slight decrease of the centroid energy and an increase of damping 
with increasing asymmetry. 

In the second case we excite protons and neutrons out of phase which excites isovector quadrupole osciUations in 
figure 2.9. The isospin energy and the isospin show nice oscillations in the quadrupole multipole moment. There is a 
slight excitation of an isovector dipole mode too. We would like to remark that due to the test particle width there 
occurs a mode coupling Is^ which can be seen in the slight excitation of isoscalar modes in the kinetic energy. The 



isospin dependence shown in figure p.7| has the same qualitative behaviour as the isoscalar mode, but is stronger in 
this case. 
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Figure 2.8: The time evolution 
of multipole moments for the 
isoscalar quadrupole excitation 
of 40^?" analogously to figure ^.2\ 
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Figure 2.9: The time evolution 
of multipole moments for the 
isovector quadrupole excitation 
of 40^?" analogously to figure p.2|. 
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Chapter 3 

Collective modes in linear response 



Now we would like to focus on small amplitude oscillations and we will develop the linear response theory. Large 
amplitude motions and effects beyond linear response are discussed in chapter ^. The simplest microscopic theory 
which provides basic experimental features and which allows to include coUisional correlations is the Fermi-gas (liquid) 
model including dissipation. We will compare the results from this linear response with the simulation of finite nuclei, 
and both will be compared with the experimental data. 

The principle of linear response is easily explained. When a system which is described by the kinetic equation 
( pT| ) is disturbed by an external potential it will react and will create a density change Sp. The connection of the 



latter to the external perturbation is called the response function, dp = xU^^'^ ■ Without mean field in (2.1) we would 
obtain the polarization 11 of the system. The mean field introduces a selfconsistency such that the relation between 
response function and polarization becomes 



n 



^=-, «7TT (3-1) 



5n ' 

where we have assumed that the mean field U is only density dependent. Therefore the potential for the response 
§u_ 

Sn 



function is V = j^. The collective mode of a system is now characterized by the condition that the denominator in 



(3.1) vanishes 

e(q,c^) = l-nq)n(q,c^)=0 (3.2) 

since then an infinitesimal small external potential can create a finite density oscillation due to the diverging response 
function. Therefore we will search for complex zeros of the dielectric (or dinuclic) function (DF), e(q, + ij) = 0, 
which provides us with the energy Q and the damping 7 of the collective mode. The strength function is then given 

by 

, , 1 Imll 
^(q''^) = ^(l_^Ren)2 + (FImn)2- ^^-^^ 

The equation of state like the isothermal compressibility can be expressed by 

1 f dn\ 1 ,. f duj 1 



t2 



dp J rp n^T q 



(3.4) 



and the fluctuations and diffusion coefficients can also be expressed via the response function. 

Two lines of theoretical improvements of the response function can be found in recent publications. The first 
one starts from TDHF equations and considers the response of nuclear matter described by a non time-reversal 
Skyrme interaction jli, 58, 5^. The other line tries to improve the response by the inclusion of coUisional correlations 



60I 6l|, Is], m and by considering multicomponent systems [|l^ 66, 6^, In Q both hues of improvements 
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have been combined into one expression. We follow this line and derive the response function from a kinetic equation 
including mean field (Skyrme) and coUisional correlations. Therefore we consider interacting matter which can be 
described by an energy functional ( mean field) E originally introduced by Skyrme ^ and the residual interaction. 
The latter we condense in a collisional integral additional to the TDHF equation. Then the response to an external 
perturbation will contain the effect of Skyrme mean field and additionally the effect of the residual interaction. 



3.1 Response function for asymmetric matter 

The system now consists of a number of different species (neutrons, protons, etc..) interacting with their own kind 
and with the others. It is important to consider the interaction between different sorts of particles if we want to include 
friction between different streams of isospin components, etc. In particular the isospin current may not be conserved 
in this way. Let us start with a set of coupled quantum kinetic equations]^ for the reduced density operator pa of a 
certain species denoted by the subscript a 



dtPa{t)=l[Pa,Ea{t) + UT]- 



dt 



,P a{t')~pY-{t') 
Tab{t - t') 



(3.6) 



where E — P"^ /2m + U denotes the kinetic as well as mean field energy operator and the external field which is assumed 
to be a nonlinear function of the density. We approximate the collision integral by a non-Markovian relaxation time 
which is derived in appendix This turned out to be necessary to reproduce damping of zero sound ||7l| , [72| . 
It accounts for the fact that during a two-particle collision a collective mode can couple to the scattering process. 
Consequently, the dynamical relaxation time represents the physical content of a hidden three-particle process and is 
equivalent to memory effects. 

Furthermore, we assume relaxation towards a local equilibrium 



/:°-(p,R,i) = 



9 



{P- 



\Pa 



£a(p-Q,(R, t))-/ia(R,t) 
Ta(R,<) 



(3.7) 



with the Fermi distribution /q. This local equilibrium will be specified here only by a small deviation of the chemical 
potential of species a ensuring density conservation |64| . The extension of the method including further conservation 
laws and specifying also the local current and the temperature can be found in p^ . 

Linearizing the kinetic equation (3^) we obtain the matrix equation for the density deviation Spb due to an external 
perturbation J/^^t 



, Sub < Sab 



6ab Cab 

n 



Ila{^^ H )ciab 



n„(c. + -)t/°,, 



(3.8) 



The matrix Cab is given by 



nc(o) 



cb 



(3.9) 



with aab = -(f^Ua the linearization of the mean field with respect to the deviation of density from equilibrium value 
caused by the external perturbation C/cxt- The partial polarization function of species a is 



^a{uj) 



dp 



fa{p+l)-faip-l) 



(27r)3 ea{p+l)-eaip- 



1 1 



■ LO - 



iO 



^ The quasiclassical Landau equation (2.1) follows from gradient expansion as 



— / -I- dptdr} - dredpf = 

at T 



(3.10) 



(3.5) 



12 



dp 



faip). 



The factor 2 in front of the integral accounts for the spin degeneracy according to pa — 2 J j2^- 

Equation (3^) represents the complete polarization of the system, because Sp ^ Ilt/oxt- It represents a matrix 



equation which is solved easily. The collective modes are given by the zeros of the determinant of the matrix on the 
left hand side of ( [3.^ ) because these are the poles of the polarization matrix. Since we took into account relaxation 
processes between all species we are able to cover current - current friction. For a two component system, e.g. neutrons 
with density p„ and protons with density pp, we write explicitly 



(1 - n^a„„)(l - nfapp) - {Dnp 



n^aprO = 



with the generalization of the Mermin polarization function p3| to a multicomponent system 



UJTa-\-i 



and the additional coupling due to asymmetry in the system 



np 



Tp Cnn iuJT'fi 



(3.11) 



(3.12) 



(3.13) 



The Dpn are given by interchanging species indices. This term does not appear for symmetric matter. Therefore we 
call this term the asymmetry coupling term further on. The result (3.11) represents the generalization of the dispersion 
relation (3.2) to a two - component system. 

It is illustrative to recover known results for symmetric nuclear matter. This is performed for the case of equal 
relaxation times Tp = t„ = r and equal deviation from the mean field ai = arm = ctpp and a2 = a-np = ctpn- Eq. 
(3.11) takes then the known Mermin form of dispersion relation [Q 



1 — (ai ± Q!2)- 



n(w- 



ujT-{-i 



n(o) 







(3.14) 



with the isovector mode ai — a2 and the isoscalar mode ai +0:2- Please note that for zero temperature the Mermin 
expression ( ^.14 ) agrees with the response function derived from taking into account multipole expansion in momentum 
of the disturbed distribution p6[ . 

We have presented a general dispersion relation for the multicomponent system including known special cases. The 
dispersion relation ( 3 . 1 1 ) is similar to the one derived recently in Q if we neglect the coUisional coupling Dnp . Also 
a more general polarization function ( 3.12| ) is presented here including collisions within a conserving approximation 
| |60| . In the following we will apply this expression to the damping of giant dipole resonances in symmetric as well as 
asymmetric nuclear matter. 



3.2 Application to nuclear matter 



Before we can solve the dispersion relation ( ^.11 ) we have to specify the wave vector characteristic of the considered 



mode. This is performed according to the Steinwedel-Jensen [[75| model. We assume that the surface remains constant 
and the density oscillation obeys a wave equation with the boundary condition that the radial velocity vanishes on 
the spherical surface with radius R — 1.13A^/^. This leads to ji{kR) = with the spherical Bessel function of order 
I = 0, 1... associated with the monopole, dipole... resonances. From this condition one obtains a connection between 
the wave vector and the radius of the nuclei or the mass number in the form 

for the i-th mode of multipolarity / where ca is the z-th zero of j;(c) = 0. 

For the monopole modes as compression modes we have no zero of first order for I = 0. For the dipole modes one 
has in first order k = 2.08/iZ which describes the giant isovector dipole resonance (IVGDR) while the ISGDR is a 
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spurious mode in first order. This would just mean an unphysical oscillation of center of mass motion. However, in 
second order k = 5.94/ R the isoscalar giant dipole resonance (ISGDR) has been observed recently ([[76[ and references 



therein). This can be considered as a density oscillation inside a sphere as we will discuss in chapter 3.6 



Besides the occurring wave vector we have also to specify the relaxation times which contain the effect of collisions. 
The dynamic relaxation time has been derived by Sommerfeld expansion |^ in appendix ^ as 



1 



1 



TabiO) 



1 



-f-V 

4 VttT/ 



1 



1 



-f-V 

2 \ttTJ 



(3.16) 



for a, b neutrons or protons respectively and dependent on whether we use a Fermi liquid or Fermi gas model. The 
Markovian relaxation time is given in terms of the cross section aat between species a and b as = ^^(JgbT'^- 
We will use as an illustrative example a slightly different mean field potential than (E^) given by |7^, pS] 



(1 + y)(Pn + Pp) - {XO + ^)Pa 



J [{Pn + Ppf ~ pV\ 



(3.17) 



with a = n,p the density of neutrons or protons respectively. The Coulomb interaction leads to an additional 
contribution for the proton mean-field 



q 



The model parameters used here reproduce the Weizsacker formula 



1 = -"^ 



02 , a^Z'^ 2 



^1/3 ^4/3 



(3.18) 



(3.19) 



with the volume energy ai = 15.68 MeV, Coulomb energy a3 = 0.717 MeV and the symmetry energy 04 = 28.1 MeV 
with the asymmetry parameter / = . 



0.0 



collision + Coulomb 

Coulomb (without collision) 
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Figure 3.1: The experimental 
centroid energies (A) and damp- 
ing rates (B) [filled symbols] of 
the giant dipole resonances vs 
mass number (data from Ref. 
|7^) together with different ap- 
proximations at T = 0. 



The energy and damping rates are now dete rmin ed by the zeros of the (Mermin) polarization function ( 3.14 ). First 
we plot the solution of the dispersion relation ( 3.11 ) for symmetric nuclear matter. In figure 34 we have plotted the 
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real and imaginary (FWHM) part of complex energy for different approximations with relaxation time (collisions) 
and with and without Coulomb mean field (3.18). In figure (A) we find that the inclusion of Coulomb effects 
reproduces the experimental shape of the centroid energies at higher mass numbers (dot-dashed line). Taking only 
collisions into account fails to describe higher mass numbers (solid line). Considering Coulomb together with collisions 
(dashed hne), the centroid energies are reduced towards the data. 

The experimental values of the damping rates are also presented versus mass number (figu re |3.l| (B)). We recognize 
that the FWHM is just twice the damping rate F = 27 which has been recently emphasized [|l5[. It has to be stressed 
that the experimental data are accessible by this FWHM. Here we have considered only collisions and have a vanishing 
Landau damping for the infinite matter model |6| . 



3.3 Comparison between simulation and linear response 

We have already stressed that the simulation of finite nuclei neglecting collisional correlations lead to a damping of 
about FLand = 2MeV for the giant dipole resonance at zero temperature. This was attributed to the surface effects of 
spherical symmetric nuclei. We consider this as the numerical value of the contribution attributed to the wall |^ at 
zero temperature. For finite temperature the will be considered an additional contribution due to thermally deformed 



shapes in section 3.5. Nevertheless it is now interesting to compare the isospin dependences of the simulation without 
collisions and the linear response including collisions but no surface effects. 



i 
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Figure 3.2: The dependence of 
isovector dipole resonance energy 
and damping on the asymmetry 
obtained from simulations and 
from linear response for A = 90. 
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In figure 3.2 we see that the damping is slightly increasing with increasing asymmetry for simulations without 



collisions while the collisional contributions cause a slight decrease of damping. This decrease is much more pronounced 
in the second order isoscalar dipole mode. The energy in turn shows the same behaviour for simulations without 
correlations and for linear response with correlations; it slightly decreases the value. 



The same qualitative feature can be recognized for the quadrupole modes in figure (3.3). However the importance 
of collisional damping is in this case much larger than the finite size effects. 



3.4 Giant resonances in excited nuclei 



We consider now the temperature depe nden ce of the GDR. We have found that the centroid energy only slightly 
decreases with increasing energy. In Fig. 3.4 the theoretical damping rates F = 27 of the IVGDR modes in ^^°Sn and 
^"^Pb are plotted as a function of temperature together with experimental values. The results of the Fermi gas model 



15 



C3-- 



> 25.0 



-€]-- 



-B- 



--a 



-QIVCQR Vlasov simulation 

ISGQR Lin. response bulli matter 

IVGQR Lin. response bulk matter 



15.0 
7.0 
6.0 
5.0 
4.0 
3.0 

2.oC3- 

1.0 

0.0 

0.00 



-€}-- 



O.IS 0.20 

asymmetry 5 



Figure 3.3: The dependence of 
isovector and isoscalar quadrupole 
resonance energy and damping on 
the asymmetry obtained from simu- 
lations and from hnear response for 
A = 90. 



(A.ll) and the Fermi Uquid model ( A.19| ) are very close and in good agreement with the data at low temperature. 
The higher temperature dependence still remains too flat compared to the experiments. The small difference between 
both models for T=0 vanishes with increasing temperature. 
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Figure 3.4: Experimental damping rates vs. tem- 
perature of IVGDR for ^^osn and ^o^Pb (^^OSn from 
Ref. Q and ^"Spb from Ref. p|) compared with 
the solution of the dispersion ( |3.14 ) relation F = 27 
for the Fermi gas (dashed lines) and the Fermi liquid 
model (solid lines). Please note that this damping 
is reduced when the momentum conservation is in- 
cluded like it is done in figure p7^. 



While the figures indicate that the linear response including collisional correlations together with the Steinwedel 
Jensen model can describe quite well the gross features of giant dipole resonances we should keep in mind that we 
have neglected here the surface effects. From the simulation results of the second chapter we have seen that there is a 
damping of about FLand = 2MeV for finite nuclei at zero temperature without collisional correlations. Simply adding 
now both contributions would overestimate the experimental damping. We would like to point out that the inclusion 
of only density conservation in the derivation of the response function so far is overestimating the width. In ||6^ we 
have shown that the inclusion of momentum conservation reduces this width by approximately 2MeV. Therefore we 
might anticipate that both the finite size effects together with the collisional contribution lead to the correct value for 
low temperatures. In figure 3.5 we see that the temperature increase is too flat compared to the experimental finding 
if we consider only collisional damping (thin lines). From simulation results of ground state we have already found a 
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• Sn, exp 
■ Pb, exp 
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Figure 3.5: The effective damping 
consisting of only coUisional damping 
compared with coUisional and surface 
damping [^^"Sn (solid lines) and -^^^Pb 
(dashed lines)] together with the exper- 
imental data (filled symbols: Sn from 
Ref. M and Pb from Ref. Q). 



contribution from finite size effects. We might now expect that the discrepancy in temperature behaviour is due to 
finite size effects. This should be explained and understood within a more simplified model. 

3.5 Nuclear surface contribution - excited nuclei 

Besides the PLand = 2MeV of the simulation in the ground state we can have additional damping due to shape 
deformation if we increase the temperature. In Ref. [p2| we have presented a method to include also scattering wit h 



the non-spherical nuclear surface. This improves the temperature dependence remarkably (thick lines) in figure 3.5 , 

The idea is to consider the surface scattering contribution to the damping determined by the additional Lyapunov 
exponent due to the surface. It has been shown that such Lyapunov exponents appear in the response function additive 
to the frequency as imaginary shift provided the Lyapunov exponent is small compared to the product of wave vector 
times Fermi velocity 

iq,R,u)^n^'{q,PfiR),u; + tX). (3.20) 

The Lyapunov exponent by itself is calculated for different deformations of the surface corresponding to different 
temperatures 

Rx{e)^Ro il + aoo + axPxicos{0))) (3.21) 

with the nuclear radius Rq — 1.13j4^/'^fm, and where A = 2 corresponds to the quadrupole and A = 3 to the octupole 
deformation . The corresponding mean deformation ^ is linked to the temperature 

J da\a\exp {—EB{a)/T) 



J da exp {—Eg [a) jT) 

where the surface dependent energy Eb(pl) is given by the Bethe-Weizsacker formula ^ 



(3.22) 



i;B(a) = -ai + ^ + -^+a4/ +05^ ^^^^ 

(3.23) 

■^For small deviations we found identical Lyapunov exponents for prolate a > and oblate a < deformations and therefore we do not 
distinguish the sign of a. 

^Please remember that in principle the Coulomb energy changes with small deformation as well according to the factor |85[ 1 — 
5(A — 1)/(2A -I- 1)^ while the surface term changes as 1 (A — 1)(A -|- 2)/2/(2A -|- Only the latter correction is considered since the 

Coulomb energy deformation would lead to corrections of around 0.3% and are neglected here. 
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with the volume energy ai = 15.68 MeV, Coulomb energy = 0.717 MeV, the symmetry energy 04 = 
the surface energy as = 18.56 MeV. For the quadrupole, S2, and octupole, S3, deformations one gets 

3(0) 5 ^^^"^ ' 

S-^ia) 5 o 



28.1 MeV and 



(3.24) 



This represents the lowest order expansion in a, however the next term gives alrea dy co rrections in fractions of percent 
for the highest deformations considered here. By this way, the statistical model ( 3.22| ) leads to a connection between 
temperature and deformation as 



(3.25) 



where the constant c is given by the coefficient of 
deformation. 



in Eq. (3.24) for the corresponding quadrupole or octupole 



Using Eq. (3.25) we can translate the Lyapunov exponent A calculated as a function of deformation into a function 
of the temperature. In figure 3.6 the contribution to the damping of IVGDR for ^^"Sn (circles) and ^°*Pb (squares) 
is presented for different shape deformations versus temperature. If we add quadrupole and octupole deformations 
we come up with a damping curve very similar to pTj . The damping starts at zero and increases rapidly with 
increasing temperature. We see that the main contribution comes from the quadrupole deformation while the octupole 
deformation is only sizeable at higher temperature. Let us note that the qualitative difference between Sn and Pb is 
reproduced by including surface scattering as well as coUisional damping. 




Figure 3.6: The scaled coUisional damp- 
ing is compared with the damping ac- 
cording to the chaotic scattering from 
the surface of quadrupole and octupole 
deformed shapes for ^^"Sn and ^"^Pb. 
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The coUisional contribution F = 27 is plotted as well (Sn: solid line, Pb: dashed line). We recognize that both 
contributions by themselves, coUisional as well as surface scattering, account almost for the same amount required by 
the experimental values, see figure 3.6. A proper relative weight between both processes is therefore necessary which 
will be introduced in the following. 

Let us note before that we have considered the surface contribution to the damping for a deformed surface by 
temperature which therefore vanishes at zero temperature. However, for a spherical surface we got already from 
simulation a zero temperature contribution of FLand = 2MeV. Consequently for the total damping caused by surface 
both contribute. 

So far we have not considered that only particles close to the surface can appreciably contribute to the surface 
scattering and to this chaotization process, while particles deep inside the nuclei are screened out of this process. 
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Consequently we consider the corresponding collision frequencies as the measure to compare surface collisions with 
inter-particle collisions. The collision frequency between particles is given by I/tq of ( 3.161 ). The collision frequency 
of particles with the deformed surface beyond a sphere, Vsmf, is given by the product of the density with the surface 
increase S{a) — S{0) — co^^ttRq according to Eq. (3.24) and with the mean velocity in radial direction Vr = 2>/%vf- 
The result is 



(3.26) 



where we have used Eq. ( 3.25| ) to replace a. We see that the frequency ( 3.26| ) is independent of the size of the nucleus 
and linearly dependent on the temperature. We use the ratio of these two frequencies to weight properly the two 
damping mechanisms, the surface coUisional, A, and the inter-particle coUisional, 7, contributions. Consequently the 
full - width - half - maximum (FWHM) reads 



Tfwhm = 2C7 + 2(l-C)(A)+27Land 



With the help of (|3T|) and ( |3^ ) the weighting factor C is given by 

C(t) 



(3.27) 



(3.28) 



One sees that for zero and high temperatures C = 1 ^tnd due to Eq. (3.27) only the coUisional contributions matter. 
Since Vsmf is linear in the temperature and 1/t depends quadratically on the temperature, the weighting factor C, 
has a minimum at temperatures around Tc = for the gas model ( |3.16 ) and the surface contributions become 

important. In the case of the IVGDR this corresponds to a temperature of T « 3.7MeV, which is the upper limit 
of currently achievable experimental temperatures. Therefore we can state that at low and high temperatures the 
coUisional d amp ing is dominant while for temperatures around Tc the surface cont ribution becomes significant. 

In figure 3^ we have compared the effective damping according to Eq. ( 3.27 ) with the experimental data. Let 
us recall that the coUisional damping value is reduced here since we include momentum conservation too for the 
Mermin polarization function. On the other side the surface damping consists now of the deformation contribution 
by temperature and by the zero temperature contribution obtained by simulation of Vlasov equation. We find a 



reasonable quantitative agreement combining all these parts. This is illustrated more in detail in Fig. 3.7 where we 
have plotted the strength func tion (^ ) for ^^OSn (LHS) and ^ospb (RHS) within the Fermi gas model ( |A.ll| ) (dashed 
lines) and Fermi liquid model ( A. 19 ) (solid lines) with the normalized data from Ref. The good overall agreement 
of the shape evolution with the experiment is again accompanied by only minor differences between the Fermi gas and 
the Fermi liquid model. 



3.6 Higher order modes - surface modes 

The existence of isoscalar giant dipole resonance (ISGDR) in nuclear matter is considered as a spurious mode in most 
text books since one associates with it a center of mass motion. The more surprising was the experimental justification 
of a giant resonance carrying the quantum numbers of a isoscalar and dipole mode p7|, u6l pq|. Consequently one 



has to consider higher harmonics as explanation of such a mode |87, g9|, |90|]. Usually this mode is associated with a 
squeezing mode analogous to a sound wave |87[ ^ ^ p8[. 

In this chapter we want to discuss the influence of surface effects on the ISGDR compression mode. We will show 
that even in the frame of the Fermi liquid model such modes can be understood. Moreover we claim that the surface 
effects are not negligible for reproducing the strength function. While we have already given a phenomenological 
approach to include surface effects we want to show now a straightforward possibility to include surface effects in the 
response function. Consequently we give first a short derivation of response function including surface effects. This 
will result in a new formula in the temperature-dependent extended Thomas-Fermi approximation |]9^. 



The starting point is again the semi-classical Vlasov equation (2.1). Instead of using now a spatial homogeneous 



equilibrium as done so far we consider now explicitly the spatial dependence. Provided we know the response to the 
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Figure 3.7: The IVGDR strength 
function in i^Ogn (LHS) and 
208pb (RHS) within the Fermi 
gas model (dashed hnes) and 
Fermi hquid model (solid lines) 
at several temperatures com- 
pared with normalized data from 
Ref. IS^. 
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external potential without self-consistent mean field, which is described by the polarization function 11 

(5n(x,i)= y"dx'n(x,x',i)C/°'=*(x',t), (3.29) 



the response including mean field, X: is given by 



X(x,x',i) =n(x,x',t)+ / dxidx2n(x,xi,^) / 7/ x(x2,x',t). 



(3.30) 



(5n(x2, t) 

Therefore we concentrate first on the calculation of the polarization function 11 and linearize the Vlasov equation 



( 2.1) according to /(p, R, t) = /o(p, IV) + 5f{p, R, t) such that the induced density variation Sp{Il, t) — J dpS f / {2TTfi)^ 
reads 



dprfx'_^_,q,, Vp/o(p, x')V,, t/-'(x', t) 



(27r;i)e 



(3.31) 



Here we have employed the Fourier transform of space and time coordinates of (2.1) to solve for 6f and inverse 
transform the momentum into the form (3.31). Comparing (3.31) with the definition of the polarization fimction 
(3.2£) we extract with one partial integration 



n(x,x',cj) = -v^- 



dpdq 



,jq(x-x') 



Vp/o(p,x') 



pq 



(3.32) 



With (3.32) and (3.3C) we have given the polarization and response functions for a finite system. 

In the following we are interested in the gradient expansion since we believe that the first order gradient terms 
will bear the information about surface effects. Therefore we change to the center of mass and difference coordinates 
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R = (xi + X2)/2, r = xi — X2 and retaining only first order gradients we get from (3.32) after Fourier transform of r 
into q 



n(R,q) = - 
dp 



dp qVp/o(p,R) 



(27r?i)3 



dp /^Vp/o(p,R) p q.Vp/o(p,R) 



pq '2 

m 

,Vp/o(^,R) 



pq 



(3.33) 



where in the last equality we have assumed radial momentum dependence of t he di stribution function /p. We recognize 
that besides the usual Lindhard polarization function as the first part of ( 3.33 ) we obtain a second part which is 
expressed by a gradient in space. The first part corresponds to the Thomas-Fermi result where we have to use 
the spatial dependence in the distribution functions and the second part represents the extended Thomas-Fermi 
approximation. So far we did not assume any special form of the distribution function. Therefore the expression 
( 3.33 ) is also valid for any high temperature polarization of finite systems. 

What remains to be shown is that the response function ( 3.30| ) does not contain additional gradients. This is easily 
confirmed by two equivalent formulations of ( 3.30| ), n~^x = 1 + Vx ^md xn~^ = 1-1- x^i which by adding yield the 
anticommutator 

[Tl-\x]+ = 2+[V,x]+- (3.34) 
This anticommutator does not contain any gradients up to second order. Therefore we have [V = 5U"^'^ /5n] 

n(R,q,c^) 



o{dl) 



(3.35) 



l-1/(R,q,c^)n(R,q,c^) 

Equation (3.35) and ( 3.33| ) give the response and polarization functions of finite systems in first order gradient ap- 
proximation. 

Now we are ready to derive approximate formulae for spherical nuclei. In this case we can assume q||R and we 
have 



n(i?,q,w) = n°(i?,q,cj) 



-dr- 



q 



1 



n°(i?,q,c^) 



(3.36) 



where Ho is the usual Lindhard polarization with spatial dependent distributions (chemical potentials, density). We 
use now further approximations. In the case of giant resonances we are in the regime of small q and Imll" ~ ui such 
that 



1 



n'^ = i-imn" 

2 



(3.37) 



since Rell*' ~ 1 — c^q^/uj'^. Within the local density approximation we know that the spatial dependence is due to the 
density p{R) — poQ{Ro — R)- Since we have for zero temperature Imll'' oc Pf{p) we evaluate 



^i^Imn" = -noS{Ro ~ i?)a„Imn° = --(5(i?o - i?)Imn" 



(3.38) 



where we assumed the density dependence carried only by the Fermi momentum. Now it is straightforward to average 
( 3.36| ) over space with the help of (3.38) 



-Ro 

n(q,c^) = / di?i?'n(i?,q,u;) 



n"(q,c.) 



-nO(q,u;) 



n" -t-n 



surf 



(3.39) 
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Consequently the surface contribution to the polarization function reads finally with (3.37) 



n-'f(q,c.) - _-l-ImnO(q,c.) (3.40) 



which is real. With (|3|), (|1^) and ( p5| ) we obtain now for the structure function 



1 imno 
- - (1 _ y(Ren" + n™>-f))2 + (yimn")2 ■ ^ • 

For small q expansion we see that the pole of the structure function becomes renormalized similar to what is known 
from the Mie mode or surface plasnion mode 

2 ^ ^ (3 42) 

1 - v^n™''f 

After establishing the structure function including surface contribution we specify the model for actual calculations. 



We choose as mean field parameterization a Skyrme force (3.17) following Vautherin and Brink [ |95[ which leads to 
the isoscalar potential 



V = U„ -U, = ^ + ^no (3.43) 
4 o 



with to = -983.4 MeV fm^, = 13106 MeV fm^, xq = 0.48 at nuclear saturation density hq = 0.16 fm'^ and the 
incompressibility of K = 318 MeV. Further we employ again the Steinwedel- Jensen model where the basic mode 
inside a sphere of radius Ro is given by a wave vector 

<7.p = (3.44) 

This would correspond to the first order isovector mode . Since this mode is spurious we have to consider the next 
higher harmonics which is 

Qtsgdr = (3.45) 
Ho 

The polarization function with this second order mode contains still contributions from the spurious mode such that 
we have to subtract this part pO[ 

n?SGDR('^) = I^%q^sadr,i0) - n°((Z,p,C^). (3.46) 

In figure |3.6| w e have plotted the experimental structure function together with different theoretical estimates 
according to ( 3.46| ) and (B.41). The inclusion of surface corrections (dashed lines) shifts the structure function towards 



the experimental values. The inclusion of collisions (dot-dashed lines), which should be wrong for isoscalar dipole 
mode due to cancellation of backscattering leads to really worse results. This supports indirectly that the mode is of 
isoscalar dipole type and surface dominated. 

3.7 Simplified model for nuclear matter situation 

The resulting wave vectors have very low values compared with the Fermi wave vector in the Steinwedel Jensen model. 



This allows us to expand the Mermin polarization function (3.12) with respect to small qvc/uj ratios and Vc the sound 
velocity. We obtain 

ma Uj{uJ + l/Ta) 

PailJ-a) = 2 A^'V3/2(2a)Ccorr (3.48) 
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Figure 3.8: The experimental 
structure function (T=0) versus 
theoretical values. The bulk 
RPA result (solid lines) is com- 
pared with the extended Thomas 
Fermi approximation (surface 
corrections, dashed lines) and 
the inclusion of collisions (dot- 
dashed lines). The latter one 
should be of less importance 
due to symmetry of isoscalar 
mode. The data suggest this 
case and support surface con- 
tributions. Circles: Normalized 
data from Ref. §91 
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where the thermal wave length is = 2ttTc' / (rnaT)^ /3/2 the standard Fermi integral and Za = e^"/^ the fugacity. 
The correction constant Ccorr = 1-22 is introduced to fit the numerical solution of the dispersion relation (3.2) with 
the approximative expansion (3.47). With the help of this expansion the dispersion relation (3.11) takes the form 



= 



.2 2 



52 



(UJ + i/Tn)Tp 



2 2 

Sp9 



c2 



{UJ + i/Tp)T„ 



(3.49) 



with the partial sound velocities c and c 



^2 



OLab- 



T /l/2(Za) /l/2(Zi.) 



(3.50) 



The dispersion relation (3.11) with the dynamical relaxation times (3.16) is a polynomial of tenth (sixth) order 
corresponding to the inclusion of memory (in)dependent relaxation times. While most of these solutions will just lead 
to parasite solutions (Re a; < 0) we will get two coupled modes, i.e. the isoscalar and isovector mode. Furthermore 
a third mode appears at extreme asymmetries and/or strong collisional coupling which we will describe in the next 
section. 



3.8 New collective mode 



Now we employ the potential ( 3.17| ) and assume different neutron and proton densities. In figure 3.9 we plot the 



isoscalar and isovector modes versus temperature for ^^Ca with a small asymmetry 5 — 0.2 as well as for ^°Ca with 
an asymmetry 5 = 0.33. The kinetic energy is linked to a temperature within the Fermi liquid model via Sommerfeld 
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expansion 



Akin 



(l + ^)5/3_^(l_^)5/3 



12 



(1 + ,5)1/3 + (l_^)l/3 



(3.51) 



This connection between temperature and excitation energy is only valid for a continuous Fermi liquid model. For 
small nuclei, the concept of temperature is questionable. Some improvement can be obtained by the definition of 



temperature via the logarithmic derivative of the density of states [ [LOO 

dp 



T- 



1 



pdE,._ 



1/2 



(3.52) 



which provides Ecx ~ j{E/A) in comparison with (3.51) for small temperatures and small nuclei. We use this 
temperature to demonstrate possible collective bulk features in an exploratory sense. Of course, the surface energy 
and shell effects cannot be neglected for realistic calculations of small nuclei. 
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Figure 3.9: The energy and damping of the IVGDR and ISGDR of *^Ca (left) and ^°Ca (right) vs excitation energy. 
Besides the isovector modes (circles) and the isoscalar modes (squares) a soft third mode appears (triangles). 

With increasing temperature the isovector and isoscalar energies decrease and vanish at a certain temperature. At 
these temperatures the damping becomes twofold because the damping of the spurious mode with negative energy 
becomes different from the physical mode. We can consider this behavior of damping as a phase transition of isospin 
demixture. At the same time a very soft mode appears due to collisional coupling which is only present in asymmetric 



matter p8|. This mode is more pronounced in the next figure 3.9 for °'Ca with an asymmetry of 6 = 0.33. 
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We see that the isovector mode does not disappear but turns over into a flat decrease with increasing temperature. 
This behaviour is coupled with the pronounced soft mode. In comparison with the more symmetric ^^Ca we see a 
different behavior of the damping where the isoscalar mode vanishes. Also the isovector mode appears unique and not 
two-fold. Now a clear transition of damping behavior for the isovector mode is recognizable which can be considered 
as a transition from zero to first sound damping. 

Besides the standard isovector and isoscalar modes we observe a build up o f a very soft mode with a centroid energy 
around 1 MeV. This mode appears due to the coUisional coupling c^^ of ( 3.50 ). When we turn off the relaxati on tim es, 
i.e. the collision integral, this mode is vanishing as well as in symmetric nuclear matter, see discussion after (3.13). It 
shows that this mode appears due to coUisional coupling of isovector and isoscalar modes. The corresponding damping 
of the crossed mode is continuously increasing with temperature. 

One may argue whether this third mode can really appear in the system. A simple consideration may convince us 
about the possible existence of such mode. Let us assume a coupled set of two type of harmonic oscillators (neutrons 
and protons) interacting between the same sort of particles with strength fc„ and fcp, respectively and between different 
sorts with Let us choose for simplicity only two neutrons and two protons. Then we obtain the coupled system 
of harmonic oscillators with frequencies = kn/m, uJ^ = kp/m and oj^^ — knp/m. The solution yields three basic 
modes in the system, i.e. = 2(w^ + w^^), 4^^^, 2{ujp + uj^p). If we neglect the different coupling between neutrons 
and protons w„p we only obtain two modes analogously to isovector and isoscalar ones. We see that the coupling 
between neutrons and protons can lead to the appearance of a third mode. 

Let us now compare the found new mode with the experimental evidence. There are some hints for a soft mode in 
^^Be [101]. The authors have observed a low lying structure at around 6 MeV excitation energy with a damping of 



around 1 MeV which has not been reproduced yet even within refined coupled channel calculations ]102{ . A standard 
explanation would give as the origin a weakly-bound single particle neutron orbital. The observed broad structure at 
6 MeV might be explained alternatively as the presented new coupled mode. The centroid energy as well as damping 
width at least seem to suggest this interpretation. 
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Chapter 4 

Nonlinear effects beyond linear response 



While most of the theoretical treatments of oscillations rely on the linear response method or RPA methods, large 
amplitude oscillations require methods beyond. In particular the question of the appearance of chaos has recently 
been investigated 1 44 , 4^, ^ . The hypothesis was established that the octupole mode is over-damped due to negative 
curved surface and consequently additional chaotic damping ^ . Here we want to discuss in which conditions 

one might observe octupole modes at least in Vlasov - simulation of giant resonances. We will consider different initial 
conditions of isoscalar giant resonances and will demonstrate that the appearance of octupole modes is dependent on 
the initial configuration which in turn demonstrates the nonlinear behavior beyond linear response. 

As a first initial condition we use the ground state distribution of coordinates while the momentum distribution 
is deformed anisotropically. We have modified the momenta in a way which corresponds to a giant octupole mode. 
The local densities and currents remain the same as in ground state. Figure 4.1 shows that at start time t — there 
is a pure giant octupole, which is damped out and a quadrupole resonance develops instead. The monopole and 
dipole amplitudes which should remain constant document the stability of simulation. In agreement with the already 
mentioned hypothesis, the octupole mode is over-damped. The figure shows the nonlinear behavior of mode coupling. 
Within the linear response the damping rate is expected to be independent of initial conditions. We choose now other 
initial conditions to show that the result is very much dependent on initial conditions. Therefore we split the nucleus 
in two parts of mass ratio 3:7 in accordance with the symmetry of the octupole oscillation and accelerate both pieces 
towards each oth er. Experimentally this might be realized as a central collision of two nuclei with corresponding 
masses. In figure 4.2 a clear quadrupole resonance appears and also a smaller octupole resonance can be seen. Both 
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Figure 4.1: The picture shows the time evo- 
lution of multipole moments \J ^^^aio{t), I = 
0, 1, 2, 3 of mass density concerning to normal- 
ization, dipole, quadrupole and octupole os- 
cillation. Excitation is due to modification in 
momentum space only. 



are damped out. Consequently there is no evidence for an over-damped octupole mode in this case. 
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In order to understand the different initializations, we split the kinetic energy into a thermal part and a collective 
part according to (p^) = ((p — (p))'^) + (p)'^- We analyze the time development of the total collective energy in the 
system 

Ecouit) ^ ^ I dve{v){vY{v,t) (4.1) 



with the mean current 



2ra 



(p)(r,i) = ^ /■ p/(p,r,t) (4.2) 

q{v) J (27rj-^ 



and the density 

e(r,t) = 



/(^/(P,M). (4.3) 
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Figure 4.2: This picture shows the time evo- 
lution of multipole moments \J^^^aiQ{t),l — 
0, 1, 2, 3 of mass density. Excitation is due to 
asymmetric splitting (0.3 to 0.7) of the nucleus 
in two pieces using a plane in coordinate space. 
Then both pieces are accelerated towards each 
other. 



In figure 4.3 the development of collective energy can be seen. There is a background of about 50 MeV due to 



and (3.39) in (4.1) one obtains 

Ecoll 



fixed correlations caused by finite width of pseudo-particles as one can see from the following estimation. Using (3.7) 

AN AN 

i=l j=l •' 
AN AN 

i=i j=i 

AN 



(4.4) 



i=l 



For simulation parameters of ^'^^Pb, po = 0.162 fm^'^, N = 75, Cr = 0.53 fm, we obtain a basic collective energy of 
54MeV. Using more test-particles would diminish this level. 

The solid line corresponding to figure 4.1 shows no initial collective energy. The exclusive initial excitation in 
momentum space without correlation in the spatial domain leads to zero initial collective energy. These correlations 
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are forming during time evolution. Of course, there is no center of mass motion, otherwise we would see just the mean 
streaming velocity. 

This situation is changed if we use the second preparation with simple momentum-space-correlations. The long 
dashed line in figure 4.3 shows initial collective correlations corresponding to figure 4.2. Since we can deposit enough 
collective energy in this case, we observe a clear octupole motion. 
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Figure 4.3: The time evolution of collective 
energy. The solid line corresponds to the exci- 
tation scheme of figure 4.1 and the long dashed 
line corresponds to the figure 4.2, respectively. 
While the first starts without collective en- 
ergy, the second one starts with maximal col- 
lective energy. 



In order to compare with |4^ we calculate the adiabaticity index 77 defined in |Q as the ratio of maximum radial 
surface velocity to the maximum particle speed. A smaller ratio denotes a more adiabatic shape changes in relation 
to the particle speed. In analogy we define such index as a ratio 



(4.5) 



of the root mean square radius speed in forward direction (opening angle i?) and the Fermi velocity. With opening 
angle 0.4 rad we obtain a maximum rj = 0.12 for figure 4.1 and rj = 0.30 for figure 4.2. This shows that we are 
essentially still in the adiabatic regime described in p^ . 

The nonlinear behavior described so far already documents that we are in a regime of large amplitude oscillations 
where linear response fails. The corresponding radius elongation in coordinate space varies about 10 %. 
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Chapter 5 



Summary 



We have investigated the giant resonance oscillations by solving the collision-less Vlasov equation as well as by linear 
response theory. 

With the help of numerical solution of the Vlasov equation we could reproduce the mass dependence of energy 
with one fixed parameter of test-particle width for both giant monopole and giant dipole modes. A multipole analysis 
was performed which has allowed us to characterize the corresponding excitation. It was shown, that the asymmetry 
influences the collective behavior. With increasing asymmetry the energy decreases while the damping increases. 
The damping due to flnite systems amount to 2MeV almost constant for all mass numbers which underestimates the 
experimental values considerably. This motivates to search for additional damping mechanism which was found to be 
due to the collisional correlations. 

The linear response can be used in a simplified liquid drop model to describe the giant resonances in asymmetric 
matter. We find that the collisional contribution as well as surface contributions are both important to reproduce the 
experimental damping of giant dipole resonances. While for ground state resonances it is sufficient to add the damping 
of finite size effects from simulation with the collisional contribution from linear response theory, for the temperature 
behavior we have needed also the deformation of the surface. The combined model between surface and collisional 
contribution, weighted properly due to their collision frequencies, is able to reproduce the experimental damping curve 
with temperature as well as the structure function. A higher order mode like the recently measured ISGDR mode has 
been described within this simple linear response model. 

We have observed that due to correlational coupling there can exist a new mode which appears besides isovector 
and isoscalar modes in asymmetric nuclear matter due to collisions. We suggest that this mode may be possible to 
observe as a soft collective excitation in asymmetric systems. The transition from zero sound damping to first sound 
damping behavior should become possible to observe for isovector modes since they do not vanish at this transition 
temperature like in symmetric matter. At a certain critical temperature the collective! isoscalar mode vanishes and 
the damping becomes two-fold. This can be considered as a phase transition of demixture. 

On the example of octupole resonances in finite nuclei the dependency on the initial configuration is demonstrated. 
The appearance of an octupole mode was shown by correlating the spatial and momentum initial excitation. It is 
possible to excite an octupole mode with sufficient collective energy deposited initially. We suggest that isoscalar 
giant octupole resonances should be possible to observe in nuclear collisions of mass ratio about 3:7 corresponding to 
octupole symmetry. This effect illustrates a mode beyond linear response. 
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Appendix A 

Derivation of non-Markovian relaxation 
time 



Here we shortly sketch the derivation of the damping of a collective mode within a Fermi gas and a Fermi liquid 
model. We will show that we get the latter one from the Fermi gas model with an additional contribution from the 
quasi-particles. 

We start with the Fermi gas where the dispersion relation between momentum and energy is given by e = /2m 
and will show later what has to be changed for a Fermi liquid where e is a solution of the quasiparticle dispersion 
relation. We will see that the contributions from the quasi-particles alone leads to the Landau formula of zero sound 



damping |103, 104] 



7 (X 



1 



V27rr 



(A.l) 



Our considerations use conveniently the Levinson equation for the reduced density matrix / which is valid at short 
time processes compared to inverse Fermi energy Ti/ef and which coUisional side has the form: 



dT 



{2.nf 1"^' '''' 



Ae(r)dr/;ij 5(Ap) (/1/2/3/4 - /1/2/3/4), 



(A.2) 



where / = ! — /, Ap = pi + P2 ^ P3 ^ Pa etc., g is the spin-isospin degeneracy and the transition probability is given 
by the scattering T-matrix. In case that the quasiparticle energies e{t) become time independent like in the Fermi gas 
model, the integral in the cos function reduces to the familiar expression Aer. We linearize this collision integral with 
respect to an external disturbance according to 



f = n + 5f 



(A.3) 



where n is the equilibrium distribution. Clearly two contributions have to be distinguished, the one from the quasi- 
particle energy and the one from occupation factors First we concentrate on the Fermi gas model where we have 
only the contribution of the occupation factors and will later add the contribution of the quasiparticle energies for 
Fermi liquid model. We obtain after Fourier transform of the time 



Here we use the abbreviation 



5h = 



n) + (5_ (Ae - VL)] {SFi + SF2 - SF3 - SF4) (n) 



2g f dp2dp3(ip4 



{2TThy 



\Tf5iAp). 



n\2TThf 



(ie2 deader 



)sm 



cos(0/2) 



(A.4) 



(A.5) 
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where the last hne appears from standard integration techniques at low temperatures. Further abbreviations are 

V 

^±(2;) — 7r(5(a;) ± i — « 7r(5(x) 

X 

SFi = -^/i(n2«3"-4 + "2"'3"4)- (A.6) 

The approximation used in the first line consists in the neglect of the off-shell contribution from memory effects. This 
is consistent with the used integration technique (A. 5). This terms would lead to divergences which has to be cut off 

H- . ^ . . . . . 

Neglecting the backscattering terms 6F2/3/4 we obtain from (A.4) a relaxation time approximations with the 
relaxation time 



1 



r(ei) Att^to 



dx2dx^dxi [5{Ax + w) + (5(Aa: — to)] (ri2n3n4 + n2'n^n4) 



(A.7) 



-A 



with uj — il/T, X — {e — fJ-)/T, A — fi/T and the time 

1 2gmT'^ 



{A.. 



To ar 

Here we have used the definition of cross section |Tp = {i-Kh^ / m)^ da / dfl and have assumed a constant cross section 
a. To calculate (A.7) one needs the standard integrals for large ratios of chemical potentials /i to temperature A — ji/T 



00 

j dx2dx^dx4n2n^nA^{l^x ± a;) = — rii(xi ± lo) [tt^ + {xi ± lj)^] 

-A 



to obtain 



1 



[27r2 + (xi +uof + (xi - ujf 



(A.9) 



(A.IO) 



Further we employ a thermal averaging in order to obtain the mean relaxation time finally]^ 



= / dxiTLirii 



-A 



r(ei) 



1 

To 



If we do not use the thermal averaging but take (A.IO) at the Fermi energy ti = tf we will obtain 



(A.ll) 



(A.12) 



We see that both results disagree with the Landau result of quasiparticle damping (A.l) by factors of 3 at different 
places 104 lQt\. We have point out that the result at fixed Fermi energy will lead to unphysical results for the 

Fermi liquid case. Therefore we consider the thermal averaged result as the physical one. 

We now turn to the Fermi liquid model and replace the free dispersion e = /2m by the quasiparticle energy Cp. 
Than the variation of the collision integral gives an additional term which comes from the time dependence of the 

^One has to use the identities valid up to o(exp[— A]) 



00 00 
j dxnn = 1; j 



dxx nn = — . 
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quasiparticle energy on the cos-term of (A_.2). We have instead of the sum of two complex conjugate exponentials in 
(A.2) an additional contribution from the linearization of the exponential 




= e 



-iAe 



t-r 

1 - i A / di[e(t) - e] 



— e 



-iAe 



t 



dt 



5f{t) 



(A.13) 



In the last line we have replaced the variation in the quasiparticle energy e(t) — e by the variation in the distribution 
function 5f due to the identity |104 



fit) - n{e) = /(i) - n(e(i)) - He) ~ 

71 H 

-n'{e{t)^e)^ — [e{t)-e] 



(A.14) 



where we assumed within the quasiparticle picture that f{t) — n{e{t)). This leads now to an additional part in the 
relaxation time which we write analogously to (A.4) 



Using again 



|) we obtain 

1 



h (5+(Ae — O) — S^{Ae + fl) 71171,2713714 — 711772713774 \ 
2 O rti77i / 



-3 ( n{xi + uj) [n"^ + {xi + uj)'^ 
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and get after thermal averaging ( A.ll| 
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Tc 



1 

To 



Taking instead of thermal averaging the value at Fermi energy (ei — £/) in ( A.16| ) we find 

1 3(7r2 + Lu^) e'^ 



Tc(e/) 



1 



2tt'^touj e'^ + 1 ' 



(A.15) 



(A.16) 



(A.17) 



(A.18) 



Here we like to point out that the Landau result (A.l) appears in (A.17) (see also in Ref. 107, |l^, 105, |TOi 
Adding now dAHl ) and ( |Al7| ) we obtain a final relaxation time for the Fermi liquid model 



1 

Tliq 



2 

To 



(A.19) 



which is the main result in this paper. It contains the typical Landau result of zero sound (A.l) except the factor 2 
in front of the frequencies. Comparing ( A.19| ) with the Fermi gas model (A.ll) we see that in the limit of vanishing 
temperature the Fermi liquid value is lower with oc 2fl^ compared to the Fermi gas cx 3ft^. Further for vanishing 
frequencies (neglect of memory effects) the Fermi liquid model leads to twice the relaxation rate than the Fermi gas 
model. The coefficient of temperature increase is than twice larger for the Fermi liquid than for the Fermi gas. If we 
consider the relaxation times at Fermi energy (no thermal averaging) (A. 12) and (A. IS) we find the same results as 
above in the limit of vanishing frequencies. For vanishing temperature only the Fermi gas ( A.12| ) coincides with the 
result of ( |A.ll ) oc 3f2^. Expression (A.18) goes to zero for T = and underlines the necessity to thermal average the 
value. 



^Here one uses [o(exp[— A])] 
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